Journal of Thermal Science Vol.26, No.1 (2017) 1-10 


DOI: 10.1007/s11630-017-0902-2 


Article ID: 1003-2169(2017)01-0001-10 


Modeling and Numerical Investigation of the Inlet Circumferential Fluctuations 


of Swept and Bowed Blades 


TANG Mingzhi, JIN Donghai , GUI Xingmin 


Aeroengine Numerical Simulation Research Center, Beihang University, Beijing, China. 


© Science Press and Institute of Engineering Thermophysics, CAS and Springer-Verlag Berlin Heidelberg 2017 


The circumferential fluctuation (CF) source terms induced by the inviscid blade force can affect the inlet distribu- 


tion of flow parameters and radial equilibrium of swept and bowed blades. However, these phenomena cannot be 


adequately described by throughflow methods based on the axisymmetric assumption. A transport model for the 


CF stresses is proposed and correlated to the distribution of circulation to reflect the effect of the inviscid blade 


force. To investigate the effect of the inlet CFs on swept and bowed blades, the model is integrated into a 


throughflow model and applied to a series of cascades with different sweep and bow angles. For swept cascades, 


the CF source terms change the distributions of incidence angles, as well as the radial equilibrium at the inlet of 


the blade passage. And the influence is enhanced as the absolute value of the sweep angle increases. For bowed 


cascades, the distributions of incidence angles are also altered. For both cases, the model can offer a good predic- 


tion of the inlet CF source terms, and prove to exert a better prediction of blade design key parameters such as 


flow angles. 


Keywords: Throughflow model; Circumferential fluctuation (CF) stresses; Transport model; Swept and bowed 


blades 


Introduction 


Fully three-dimensional (3D) Navier-Stokes (N-S) si- 
mulations have been widely used in the design process of 
turbomachines nowadays. However, throughflow method 
remains an important tool for designers [1] in the early 
design process, as it can rapidly provide the predictions 
of performance and flow fields [2, 3]. 

As a dimension reduction method, the circumferential 
average method was first proposed by Smith [2]. Later, 
throughflow models based on circumferentially averaged 
Euler [4] or N-S [5-8] equations were widely investigated 
and applied to predict the performance of compressors, 
and these models were proved to be well suited for ob- 
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taining transonic compressor flow fields [9]. In these 
throughflow models, the axisymmetric assumption is 
adopted. With this assumption, the throughflow solvers 
may fail to reproduce some influences induced by 3D 
blade design, such as the change in the radial distribution 
of blade loading. 

CF stresses, related to the blade-to-blade spatial fluc- 
tuations, arise in the circumferential averaging of the 3D 
N-S equations, due to the non-linearity of the equations. 
These terms are usually neglected in most throughflow 
models, leading to the simple axisymmetric form of the 
equations. The reason to neglect thenon-axisymmetric 
character is that the magnitude of the CF stresses was 
found to be small, compared with other terms in the 
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2 
Nomenclature 
X, r axial and radial coordinate (m) 
o circumferential coordinate (rad) 
U vector of conservative variables 
F, G vectors of convective fluxes 
F,, G, vectors of viscous fluxes 
S vectors of source terms 
Fp inviscid blade force 
Fr viscous blade force 
P circumferential fluctuation source terms 
p density (kg/m”) 
w relative velocity (m/s) 
v absolute velocity (m/s) 
p pressure (Pa) 
E total internal energy (J) 


governing equations [2]. Another reason is that the de- 
signers do not have a priori knowledge of their magni- 
tudes However, the CF stresses might be more important 
than the viscous terms, and their importance grew with 
the loading [10, 11]. With the CF terms included, a major 
enhancement in the prediction capability of the through- 
flow analysis is obtained. Besides, the CF stresses might 
affect the spanwise and axial distribution of flow para- 
meters, such as velocity, flow angle, pressure and loss [8, 
12, 13]. These stresses can reveal the effect of corner stall 
and radial mixing. Moreover, the CF stresses will also 
alter the distributions of flow parameters at the inlet of 
the blade rows. For swept or bowed blades, the radial 
equilibrium would change, resulting in a stronger redi- 
stribution of flow parameters and reflecting the effect 
induced by swept or bowed blades [14]. 

The CF stresses are mainly induced by the tip leakage 
vortex, passage shock, shock-boundary layer interaction 
and hub corner stall [10], as well as the inviscid blade 
force [14]. However, the difficulty still exists in modeling 
these physical phenomena to obtain the CF stresses. The 
CF stresses exhibited a harmonic behavior [15] and were 
modeled by a nonlinear harmonic method [16]. However, 
this approach was difficult to apply to the blade-to-blade 
stream surface. 

In this paper, a stress transport model is formulated to 
calculate the CF stresses. This model aims at providing a 
quantitative approach to account for the nonaxisymmetry 
of the flow. The effect of CF stresses on swept and 
bowed blades are also investigated in the paper. 


Governing equations 


Circumferentially averaged N-S equations 
The throughflow model is based on the circumferen- 
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Le total enthalpy (J) 

w rotating speed (rad/s) 

T viscous stress tensor 

q heat flux 

k ratio of specific heat 

Superscripts 

i circumferential fluctuation 

ü circumferential fluctuation (density-weighted) 
circumferentially averaged parameter 

_ circumferentially averaged parameter 
(density-weighted) 

Subscripts 

X, r, u axial, radial and circumferential components 

S suction surface 

p pressure surface 


tially averaged N-S equations. The circumferentially av- 
eraging operation of a flow variable q is defined as 


= 1 po 
q (x, r,t) = ean q(x, r, p,tdo 


Ap 

where Ag is the azimuthal width of the blade passage. 
For compressible flow, the density-weighted circumfe- 
rentially averaging operator is defined as 


ct (2) 
p 

Further, the 3D flow parameter is decomposed into the 
sum of a circumferential average and a spatial fluctuation 
which arises from the flow non-uniformity in the cir- 


cumferential direction, which can be expressed as 


(1) 


q(x,r,p,t)=q(x,r,t)+q"(x,r,p,t) (3) 
Due to the non-linearity of the equations, the fluctua- 
tions appear in the form of the circumferentially averaged 
stresses, i.e. the so-called CF stresses. These stresses can 
be considered as correction terms in the momentum and 
energy equations which reintroduce the mean effects of 
the flow nonaxisymmetry into the axisymmetric equa- 
tions. The non-linear terms can be expressed as 


nn 


PNI2 = PUNI + PNA (4) 
Applying the circumferentially averaging operator to 
the N-S equations, the circumferentially averaged N-S 
equations in the relative frame of reference are formu- 
lated with the CF stresses being grouped as source terms, 
i.e. the so-called CF source terms P, 
oU 10 oa 1 ô = ta 
= +55. 1or(F-F) |+5-— [or (G-G,) | © 
= S+ F, + F;+P 


where 
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In this set of equations, Fg is the inviscid blade force 
to model the effect of pressure gradient on flow turning, 
while Fr is the viscous blade force to model the blade sur- 
face friction effect on the passage flow, and b is the bloc- 
kage factor reflecting the blade geometric blockage effect. 

To investigate the effect of inlet CFs on swept blades, 
the radial momentum equation for inviscid flow can be 
induced to 


= = —?2 
1 lbp 
— ) Yu Ym coso 
bo or r i 
RGP CENT _W CENT _M (6) 
ee Wry i Fp, 
—w,, sino ++ 
dm, p P 
DWMDM 


In equation (6), RGP is the radial gradient of pressure. 
CENT_W is the centrifugal force induced by the absolute 
circumferential velocity. CENT_M is the centrifugal 
force due to the curvature of the meridional streamline. 
DWMDM results from the acceleration of the meridional 
velocity along the meridional streamline. Before the 
leading edge of the blade passage, the radial component 


of the inviscid blade force Fg, is 0, while Pr (actually the 
Pr in equation (5) divided by the averaged density) still 
exists to affect the radial equilibrium, acting as an exten- 
sion of the inviscid blade force. 


Transport equations for CF stresses 


To model the CF stresses wiw; in equation (5), the 


transport equations for the CF stresses are derived by 
analogy with the Reynolds stress-transport model, and 
presented here as: 


aloramnivt) oler) 
brox, 7 brox;, 
Sw’ 
CONV TRAN 
-== ôw, -== ôw, 
= ww; Jy n" w” 1 (7) 
a OX, ia oa 
MM. xr—q_ eur 
GEN 


The transport equations consist of six equations, but 
the transport equations for the velocity-enthalpy correla- 
tion terms in the energy equation are not modeled, since 
they are less important [17]. 

In equation (7), term CONV is the convection term, 
representing the convection of CF stresses. Term TRAN 
is the transport term. Term GEN is the generation term, 
representing the coupling of the transport equation and 
the axisymmetric flow. Term PRE is the pressure-velocity 
correlation term, representing the work done by the 
pressure fluctuation. Term CENT is the velocity-source 
correlation term, representing the work done by the 
source terms of the mean flow. 

Different from the Reynolds stress-transport models, 
the viscous terms are neglected in equation (7). The rea- 
son is that the focus of the present paper is to model the 
CF at the inlet of the blade passage, where the loss is 
small. Furthermore, CFs of the first half of the blade 
chord mainly comes from inviscid effects [11]. 

The PRE term is found to be the dominant term at the 
right-hand side of equation (7). So only the PRE term is 
modeled to close the transport model, while the other 
terms are neglected. 


Modeling the pressure-velocity correlation 


The method to model the pressure-velocity correlation 
is mainly based on the inviscid analysis. Since the pres- 
sure potential field dominates the circumferential non- 
uniformity, considering the Helmholtz decomposition, 


v=Vo+Vxy (8) 
For modeling the circumferential non-uniformity at 


the inlet of the blade passage, the flow can be assumed to 
be irrotational in the absolute frame of reference. That is, 


> 


Vxv=0 (9) 
This reduction of the problem omits the influence of 
the secondary flow which is another source of the CF 
stresses. Nevertheless, it does not introduce a lack of 
generality, since the core-flow are far from being neglig- 
ible [8, 11, 15]. 
Besides, a further assumption is that the flow variables 
p, p and w; are able to be expressed through a Taylor se- 
ries expansion in the circumferential direction. In the 
present investigation, the flow parameters are expressed 
through Taylor series expansion with partial derivatives 
up to 1“ order, i.e. 


— 0 
p=pt+t“(e-) 
OP 


— ô 
p=ptL£(o-%) (10) 
Op 
= ow, 
Ye oor (p-p) 
Op 
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which means that these variables are assumed to be 
linearly distributed in the circumferential direction. This 
assumption may fail if a strong shock wave exists in the 
flow field. However, this situation is not considered in 
the present investigation. Besides, higher order Taylor 
series expansion may bring higher precision, but the im- 
provement in precision is not investigated either. 

From equation (10), the PRE term can be derived as 


, op ôw, 6 — 6 
wr? = (9-9)-— p+ (w-e)] (D 
Ox; OY Ox ; Op 


j 
If we further assume that @ locates at the middle of 


Í 
the blade passage, and a ; a2 
op ox Op 


o, then the PRE term can be derived as 


are independent of 


pP Sh OP AG OP (9, gy) 
Ox; Op 0x,0g 12 Op Ox; (12) 
Ow, Op O 

Pict eck =) 


Op OY Ox; 


According to equation (2), the relationship between @ 
and g is 


= TD (13) 


As long as we obtain the circumferential derivatives of 
pressure and velocity, the pressure-velocity correlation is 
closed. Moreover, to formulate the correlation with the 
distribution of blade loading, which reflects the inviscid 
blade force, the formula in Wu’s [18, 19] hub-to-tip 
stream surface theory is adopted, i.e. 

Op _ DV, 7) 
0p * Dx 

So the circumferential derivative of pressure is corre- 
lated to the distribution of the circulation, which is a key 
parameter in the design process. 

To get the circumferential derivatives of the three ve- 
locity components, the continuity equation and the invis- 
cid form of energy equation are adopted. Combined with 
equation (9), the circumferential derivatives of the three 
velocity components are as follows, 


(14) 


ow, _ Ow, 

0p Ox 

Ow. Ow,r Ov, r 

— = —*- + 20r = —— 

0p or or 

é a(rpw,) aew) 
Wu yeep (Ed PW h 

OP F Ox 


Ow Ow 
—pw, (Ww, —*+ w, “)\/Le(w;, —kRT)] 
Op Op 
Considering equations (14) and (15), the major 
sources of circumferential non-uniformity are the partial 
derivatives of the circulation. 
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Results and discussion 


A series of swept and bowed cascades are selected to 
validate the model and investigate the effect of the inlet 
CFs. The swept and bowed cascades have the same ele- 
ment profile. The thickness distribution of the CDA pro- 
file is the same as SKG3.6 [20]. The basic parameters are 
listed in table 1. 

For all cases, the nominal incidence angles are zero, 
which means the inlet flow angle is equal to the blade 
geometric inlet angle. 


Table 1 Basic parameters of cascade 


Parameters Values 
Chord 100 mm 
Blade geometric inlet angle fb), 40° 
Aspect ratio of cascades 6.0 
Solidity 2.15 
Camber angle 40° 


Both 3D and throughflow numerical simulations are 
carried out. The 3D numerical solver is FINE Turbo by 
NUMECA International. The number of grid nodes of 
each cascade is 1623296. Fig. 1 shows the grids for 3D 
simulations. Two sets of Euler throughflow simulations 
denoted as CAM, representing the original throughflow 
model, and CAM+MODEL, representing the through- 
flow model including the transport model, respectively in 
the following sections, are also carried out. For CAM+ 
MODEL, the CF source terms are taken into account 
before the leading edge of the blade row, while they are 
set to zero at other portion. Total pressure, total tempera- 
ture and the direction of the flow are imposed at the inlet 
of the computation domain. The mass flow is imposed as 
the outlet boundary condition to ensure the same inlet 
Mach number. For 3D simulations, the adiabatic non- 
slipping conditions are applied to the solid boundaries. 


— 


Fig. 1 3D computational grids of the cascades 


In the following sections, the data in the figures is ex- 
tracted at about 4% of the chord before the leading edge. 


Cascades of different sweep angles 


To investigate the effect of the inlet CFs on swept 
blades, a series of simplified cascades with different sweep 
angles are employed to be compared and analyzed. In 
this paper, we use the definition of sweep defined by 
Ramakrishna and Govardhan [21], as shown in Fig. 2. 
The sweep angles are +5°, +10°, +15°, +20°, +25°, where 


the negative and positive values represent backward (BS) 
and forward sweep (FS) respectively. Fig. 3 illustrates the 
meridional view of the forward and backward cases. 


Sanaa 
Axial Direction 


Chord Line 


Fig. 2 Definition of sweep [22] 


Sweep angle (degree) 
FS20 (+20) 
FS10 (+10) 
BS10 (-10) 
BS20 (-20) 


Flow direction 
—_ > 


Fig.3 Sketch of swept cases in meridional view 


For 3D numerical simulations, viscous calculations are 
carried out. Therefore flows near the endwalls are af- 
fected by the boundary layer. In this paper, only the in- 
viscid phenomena are discussed, so the parameters near 
the endwalls are not included. 

Fig. 4 shows the relationship between the incidence 
angle and the sweep angle. At 20% height of the cascades, 
the incidence angle increases with the increase of the 
sweep angle. While it reveals an opposite trend at 80% 
height. Besides, compared with the 3D results, CAM+ 
MODEL gives a better prediction of incidence angles 
than CAM. And the difference ranges from 0.2° to 0.5°. 
This remark shows that the throughflow model without 
the CFs is not sufficient to have an accurate description 
of the flow turning. 

The improvement in the prediction of incidence angles 
mainly results from the inclusion of Pu, Pr and Px, which 
affects the distribution of meridional and circumferential 
velocity. Fig. 5 presents the relation between the sweep 
angle and Pu before the leading edge of the cascade. 
From the definition, the inviscid blade force is 0 before 
the leading edge. However, Pu still exists and will work 
as a source term to change w, along the streamline direc- 
tion. Pr and Px will affect the distribution of pressure and 
velocity in the radial and axial direction respectively. As 
stated in [14], the CF source terms are induced by the 
inviscid blade force, and work as the extension of the 
inviscid blade force to the upstream of the blade passage. 


—#— 3D 


—e— CAM 


—4— CAM+MODEL 


0.8 


o 
P 


incidence angle 
2 
oO 


ò 
IN 


ò 
for) 


-30 -20 -10 0 10 20 30 
sweep angle 


(a) 20% 


—s— 3D 
—e— CAM 
—4— CAM+MODEL| 


incidence angle 


sweep angle 


(b) 50% 


incidence angle 
I 


-30 -20 -10 0 10 20 30 
sweep angle 


(c) 80% 


Fig. 4 Relation between the sweep angle and the incidence 
angle at different spanwise position 


(2s Pp) 


rA 
at the present working conditions of these cascades. 
Therefore, Pu is also negative, representing the extension 
of pressure difference on the blade surfaces at the leading 
edge of the cascades. As observed in Fig. 5, the transport 
model can give a good prediction of Pu. The magnitude 
of CAM+MODEL is the same as 3D. Moreover, CAM+ 
MODEL gives a similar trend in the changes of Pu with 
sweep angles to the 3D results. The absolute value of Pu 
of both numerical simulations increases as the absolute 
value of the sweep angle increases, revealing stronger 
CFs induced by blades of higher sweep angles. As a re- 
sult, the modification of v„r at the inlet of the blade pas- 


For stators, F’z, equals , which is negative 
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Fig. 5 Relation between the sweep angle and the inlet Pu at 
different spanwise position 


sage by the CFs is enhanced. And the model can provide 
an accurate description of this modification. Moreover, 
the change differs at different spanwise positions for a 
swept blade, or for blades of different sweep angles, and 
the transport model can give a good prediction of the 
change in v,r. 

To further exploit the improvement brought by the 
transport model, the prediction of the radial equilibrium 
at the inlet of the blade passage is analyzed. The span- 
wise distributions of the terms in equation (6) are ex- 
tracted, as shown in Figs. 6-10. From these figures, 
CENT_M and Pr are the major terms to balance RGP, 
while CENT_W and DWMDM are about 10° to 10° 
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order of RGP. Hence, it is necessary to model Pr to get a 
more accurate prediction of the radial gradient of pres- 
sure. Furthermore, for 3D, from BS10 to BS20, the value 
of RGP increases from about 5000 m/s* to about 12000 
m/s? (Fig.6), and the value of Pr increases from about 
1000 m/s? to about 7000 m/s? (Fig. 10), which makes 
most of the contribution to the change in RGP, while 
CENT_M increases from about 3000 m/s* to 5000 m/s’. 
For forward sweep cases, similar trends can be observed. 
That is, as the absolute value of the sweep angle gets 
larger, the absolute values of RGP and Pr increase faster 
than CENT_M. Moreover, Pr gets more prominent as the 
sweep angle increases. 

From BS10 to BS20, the discrepancy of RGP between 
3D and CAM increases from about 3000 m/s* to about 
6000 m/s’, while CAM+MODEL represents a better de- 
scription of RGP than CAM, although discrepancies still 


exist between CAM+MODEL and 3D, as shown in Fig. 6. 


Besides, as the absolute value of the sweep angle in- 
creases, a more accurate prediction of RGP by CAM+ 
MODEL is obtained. Combined with Figs. 7-10, the im- 
provement mainly comes from Pr. Thus the inclusion of 
Pr by the throughflow model may provide a more accu- 
rate prediction with the increase of the absolute value of 
the sweep angle. 
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Fig. 6 Spanwise distribution of the inlet RGP with different 
Sweep angles 
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Fig. 7 Spanwise distribution of the inlet CENT_W with dif- 
ferent sweep angles 
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From Fig. 8, CAM+MODEL can also give a good de- 
scription of CENT_M. Since the difference in the meri- 
dional velocity between CAM, CAM+MODEL and 3D is 
small, CAM+MODEL can better depict the streamline 
curvature at the inlet of the blade passage, which is mod- 
ified by the new radial equilibrium. 
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Fig. 9 Spanwise distribution of the inlet DWMDM with dif- 
ferent sweep angles 


0.9 


0.8 


0.7 


0.6 


 idddddd 


--O-- FS20_3D 


0.5 


—E— BS20_CAM+MODEL| 


spanwise 


0.4 --C-- FS20_CAM+MODEL 


DD PPPPPPPPPPPPPPPPP PPP Pith 


0.3 


DD DDD DDDDDDPDDDDDDDDPPPPDDDDDSB Bp 


0.2 


DodoococooocooRooo0oRooO cococcccoorin 
Do0 S60 cc0c0coc000ccc000000e00000000g0 


0.1 
-6000 -4000 -2000 0 2000 4000 6000 8000 


Pr 


Fig. 10 Spanwise distribution of the inlet Pr with different 
sweep angles 


For FS10 and BS10, the value of Pr predicted by the 
transport model is about 0.08 of the 3D results (Fig. 10), 
and this ratio increases to about 0.2 for FS20 and BS20. 


Therefore, the transport model still needs to be modified. 
However, the model can already give a good prediction 
of the radial equilibrium at the inlet of the blade passage. 
Besides, better prediction of Pr can be obtained by the 
transport model with the increase of the absolute value of 
the sweep angle. 


Cascades of different bow angles 


To further investigate the effect of inlet CFs, a series 
of simplified bowed cascades with different bow angles 
are employed to be analyzed and compared. In this paper, 
we use the definition of bow defined by Wang Zhongqi 
[23]. The bow angles are +5°, +10°, +15°. An axial view 
of the bowed cases is presented in Fig. 11. 

The relation between the incidence angle and the bow 
angle is analyzed, as shown in Fig. 12. Similar to the 
Swept cascade cases, CAM+MODEL also provides a 
better prediction of the incidence angles than CAM. Be- 
sides, the incidence angles vary as the bow angles change, 
and CAM+MODEL can better depict the trend at 50% 
and 80% height of the cascades. 

Fig. 13 is the relation between the inlet Pu and the 
bow angle. As aforementioned, for bowed cascades, Pu is 
also induced by the inviscid blade force. Thus its value is 
negative. For bowed cascades, the transport model can 
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give a good prediction of Pu as well. 

From Figs. 12 and 13, for negative bowed cascades, 
the incidence angle decreases as the bow angle increases, 
reflecting an unloading process, thus inducing the de- 
crease of the absolute value of Pu. For positive bowed 
cascades, the process is on the contrary. And CAM+ 
MODEL can depict this trend as well. 

The radial equilibrium of the bowed cascades case is 
also altered. However, the change of terms in equation (6) 
is small. Therefore, they are not presented in the present 


paper. 
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Axial view 


Fig. 11 Sketch of bowed cases in axial view 
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Fig. 12 Relation between the bow angle and the incidence angle at different spanwise position 
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Fig. 13 Relation between the bow angle and the inlet Pu at different spanwise position 
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Conclusions 


The CF stresses are always neglected in throughflow 
models, which may fail to describe some flow phenome- 
na induced by 3D blade design. To conquer the short- 
comings, a transport model for the CF stresses is formu- 
lated and integrated into a throughflow model. The 
transport model is based on potential flow analysis and is 
correlated with the distribution of circulation, which re- 
flects the inviscid blade force. 

To investigate the effect of the inlet CFs on swept and 
bowed blades, such as the change in flow angles and the 
radial equilibrium, a series of swept and bowed cascades 
are employed. Some conclusions are summarized as fol- 
lows: 

1. The CF source terms can offer a better description 
of the radial equilibrium at the inlet of the blade passage 
induced by swept and bowed cascades. 

2. The CF source terms have a greater influence on the 
inlet radial equilibrium of swept cascades than bowed 
cascades. Besides, as the absolute value of the sweep 
angle increases, the influence of CF source terms gets 
stronger. And the accuracy of the transport model im- 
proves meanwhile. 

3. Compared with 3D results, throughflow models 
with the transport model included can provide a better 
prediction of the flow fields of swept and bowed cas- 
cades. Moreover, the transport model can provide a good 
prediction of the inlet CF source terms, i.e. a good quan- 
titative approach to account for the nonaxisymmetry of 
the flow. 
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